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We discuss statistical and physical properties of cosmic microwave background polarization, 
both in Fourier and in real space. The latter allows for a more intuitive understanding of 
some of the geometric signatures. We present expressions that relate electric and magnetic 
type of polarization to the measured Stokes parameters in real space and discuss how can be 
constructed locally. We discuss necessary conditions to create polarization and present maps 
and correlation functions for some typical models. 



1 Introduction 

Cosmic microwave background (CMB) anisotropies offer one of the best probes of early universe, 
which could potentially lead to a precise determination of a large number of cosmological param- 
etersi'i'i. The main advantage of CMB versus more local probes of large-scale structure is that 
the fluctuations were created at an epoch when the universe was still in a linear regime. While 
this fact has long been emphasized for temperature anisotropies (T), the same holds also for 
polarization in CMB and as such it offers the same advantages as the temperature anisotropies 
in the determination of cosmological parameters. The main limitation of polarization is that it 
is predicted to be small: theoretical calculations show that CMB will be polarized at 5-10% level 
on small angular scales and much less than that on large angular scales. Future CMB missions 
(MAP, Planck) will have sufficient sensitivity that even such low signals will be measurable. 
This will allow one to exploit the wealth of information present in the polarization. 

Recent work has emphasized the rich geometrical structure present in polarization Bill. In 
particular, linear polarization has been decomposed into electric {E) and magnetic {B) types, 
which transform as scalars and pseudoscalars, respectively. With polarization there are three 
additional power spectra that can be measured, in addition to E and B autocorrelation there is 
also E and T cross-correlation. Each of these can provide unique information about our universe. 



Most of this work has developed the formaUsm by using multipole expansion on a sphere or on 
a plane. Here we will develop some of the properties of polarization fields E and B directly in 
real space, which allows for a more intuitive understanding of their geometrical properties. 



2 Statistical description of polarization 

The CMB radiation field is characterized by a 2 x 2 intensity tensor lij. The intensity tensor 
is a function of direction on the sky n and two directions perpendicular to n that are used to 
define its components (61,62). The Stokes parameters Q and U are defined as Q = (/n — 122)/^ 
and U = /12/2, while the temperature anisotropy is given by T = {In + /22)/4 (the factor of 4 
relates fluctuations in the intensity to those in the temperature) . In color plates we represent the 
polarization using "vectors" with magnitude P = ^/Q"^ + U'^ that form an angle a = \ arctan(^) 
with 61. These are not real vectors since they map into themselves after 180° rotation, which 
is why we do not assign them a direction. In principle the fourth Stokes parameter V that 
describes circular polarization would also be needed, but in cosmology it can be ignored because 
it cannot be generated through Thomson scattering. While the temperature is invariant under 
a right handed rotation in the plane perpendicular to direction n, Q and U transform under 
rotation by an angle as 

Q' = Qcos2il; + Usm2ip 

U' = -Qsin2V' + ?7cos2^ (1) 

where e'^ = cos V' ei + sinV' §2 and §2 = — sin^ ei + cosV' §2. The Stokes parameters are not 
invariant under rotations in the plane perpendicular to n. For this reason it is more convenient 
to work with scalar and pseudoscalar polarization fields E(n) and B{n), respectively, which 
are invariant under rotations. In the small scale limit n is close to z and we can parametrize 
the direction in the sky with two-dimensional angle relative to a fixed coordinate system 
perpendicular to z. The two rotationally invariant fields can be written in terms of the Stokes 
parameters as 

E{1) = J cPe [Q{e) cos{2(j)i) + U{0) sm{2(Pi)] e''^ ^ 

B{1) = J (fe [U{e)cos{2(t)i) -Q{e)sm{2(Pi)] e-'''-^, (2) 

Here E{1) and B{1) are components of the two scalar fields in Fourier space. To obtain them in 
real space we can perform a Fourier transform 

E{e) = (27:)-^ J (fl e'^^ E(l) 

B{e) = (27r)-2 j dH e'^^ B{1). (3) 

These two quantities describe completely the polarization field. They can be expressed directly 
in terms of real space quantities Q{e) and U{e) as 

E{e) = -J d^e' uiO) [Q{e') cos(2^) - U{e') sin(2^)] 

= - j d^e' u{e) Qr{e') 

B{e) = -J d^e' Lo{e) [U{e') cos(2^) + Q{e') sin(2^)] 

= - j d^e' u{e)Ur{e'). (4) 



The variables {9, cj)) are the polar coordinates of the vector 6 — 6'. In equation Qj. and Ur are 
the Stokes parameters in the polar coordinate system centered at 0. For example, if 6 is zero, 
Qr = cos 2^' Q{e') - sin 20' U{e') and Ur = cos 20' U{e') + sin 20' Q{e'). The window can be 
shown to be u{e) = I/tt^^ ^ q)^ = (0 = 0). 

From equation Q we can understand how to obtain two rotationally invariant quantities 
out of Q and U directly in real space: to obtain E{0) and B{0) we radially integrate over 
circles centered at the values of Qr and Ur, respectively. Each circle is weighted by ijj{9). By 
construction these two quantities are rotationally invariant: the Stokes parameters Qr and Ur 
do not depend on the coordinate system, as they are defined relative to the — 9' vector. The 
weight function oj is also rotationally invariant. We are giving the same weight to all the points 
in each circle and we are using the Stokes parameters defined in their natural coordinate system, 
defined parallel and perpendicular to the line joining the two points (denoted with r in equation 

The variable B is clearly a pseudoscalar because it is the average of Ur and Ur changes sign 
under parity. Note that the window extends out to infinity and so the quantities E and B are 
nonlocal. This is not the only possible choice: one can construct finite extent versions of E and 
B from Q and U using a compensated window such that / OdO g{9) = 0, e.g. E = J Eg{9)d?6 
and similarly for B. The relation between these averaged E and B and Q, U is still given by 
equation ^ using 

u;{e) = -g{e) + ^jye'e'g{e'). (5) 

We can therefore construct pure E or B quantities which are obtained from an integral over 
a finite region and test the geometrical structure of polarization without measuring the whole 
sky, as long as the measured field is contiguous. Note that for these real space quantities there 
is no need to worry about mode coupling because of incomplete sky coverage, which would for 
example mix some E into B in Fourier space. 

Let us further explore properties of E and B type polarizations in real space. For simplicity 
we will continue to use small scale limit. Plate 1 shows the E polarization field and the polariza- 
tion vectors for a typical (standard CDM) model. There is only E type polarization associated 
with this model, because this is the only pattern that is produced by density perturbations 
and we did not include a stochastic background of gravity waves. We discuss this point further 
below. In figure |l| and plate 1 we can see that hot spots of the E map correspond to points 
with tangential polarization pattern (negative Qr)- Radial polarization pattern is found around 
the cold spots of E. From equation § we can read directly the relations between the pattern 
of polarization and the sign of E and B. Note that there is a formal similarity between weak 
lensing and polarization: E plays a role of convergence k, while Q and U can be viewed as the 
two shear components 71 and 72. Since weak lensing is only important for scalar perturbations 
B mode is not excited, which provides a consistency check. The sign convention for E agrees 
with weak lensing 0: positive mass density (positive k) generates a tangential pattern of shear, 
so likewise positive E gives rise to a tangential polarization pattern. 

To obtain B type polarization we can rotate all polarization "vectors" by 45°. Hot and 
cold spots of the B field correspond to places where the polarization vectors circulate around 
in opposite directions (figure |l| and plate 2). It is clear from this figure that such polarization 
pattern is not invariant under refiections (parity transformation). This is the main distinction 
between E and B type of polarization: under parity operation E transforms as a scalar and B 
as a pseudoscalar. 

We can now address the polarization pattern induced by different types of perturbations. 
In particular, we want to show why scalar perturbations cannot induce B type of polarization. 
Without loss of generality we may consider only one Fourier mode of scalar perturbations at a 
time, since an arbitrary field can always be expanded in a linear superposition of Fourier modes. 
Let us assume we are observing in direction n, while the mode is oriented in direction k. In the 





Figure 1: Polarization patterns that lead to positive and negative values of the E and B fields. The Stokes 
parameters are measured in the polar coordinate system centered at the cross. All four patterns are invariant 
under rotation but the two patterns that generate B are not invariant under reflections. 



plane perpendicular to n we can define a local coordinate system with axis aligned parallel and 
perpendicular to the plane that contains both k and n. Scalar modes have a rotation symmetry 
around k, which implies that [/ = in this frame (figure ^). Let us now look at the polarization 
pattern in the circle around direction n in the plane perpendicular to it. Polarization is directed 
either in the direction towards k or perpendicular to this direction, since only Q is generated 
in the frame determined by that direction (for simplicity only radial case is depicted in figure 
|2[). Polarization is invariant under the refiection across the axis determined by the k and n 
line, because polarization amplitude only depends on the angle between k and n. Therefore any 
integration around this circle will produce only E and not B. Since both E and B are obtained 
by radial averaging of contributions from individual circles and B vanishes for each for each of 
them it follows that i? = in general. By construction this is coordinate independent, hence 
one cannot produce B polarization with scalar perturbations in general. 

For tensor perturbations this argument no longer applies, as shown in figure H. Here the 
amplitude of perturbation depends also on the orientation of the stretching and squeezing of 
space in the direction perpendicular to the direction of plane propagation (figure ||). There is 
no parity symmetry and both E and B will be generated in general. This is also true for vector 
modes. For these modes velocity is perpendicular to the direction of wave propagation. Its 
direction in the plane perpendicular to k breaks the symmetry and generates B modes. 

3 Generation of Polarization 

The standard picture of big bang model is that the universe was hot in the beginning and was 
expanding and cooling thereafter. At high temperatures electrons and protons were free in the 
universe, because photons were sufficiently energetic to prevent any recombination to survive. 
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Figure 2: A single mode of density perturbations has symmetries under rotation around the k axis and reflection 
around any plane containing it. Only Q Stokes parameter can be present in this reference frame, implying no 
B polarization is created (see the text). Gravitational waves do not have these symmetries as illustrated by the 
deformation suffered by a ring of test particles as wave propagating in direction k passes by. There is a particular 
direction associated with the direction of the stretching and there is no rotational or parity symmetry. The same 
is true for vector modes, which induce velocities perpendicular to the wave direction. 



Once the energy of most of the photons dropped below 13.6eV ionization was no longer possible 
and at that time (at energies around 0.3eV) the universe recombined and became transparent for 
photons. The recombination was very rapid. Before recombination the density of free electrons 
was very high, so that the photon mean free path was very short. Together with electrons 
they formed a single fluid, with pressure provided by the photons and inertia by the baryons. 
This fluid supported the analog of acoustic oscillations where both the density and velocity were 
oscillating functions of time&0. The density was proportional to cos(cs/cr + (^), while the velocity 
to sm{cskT + (5), where Cg is the sound speed and the phase 5 depends on initial conditions and 
is in general time dependent. 

To generate polarization we need Thomson scattering between photons and electrons. Elec- 
tric field incident on the electron causes this electron to oscillate in the direction of the field 
and electron radiates according to the dipole emission formula. Looking for example perpendic- 
ular to the incoming light direction only one polarization of incident light will be seen to cause 
electron to oscillate, the one perpendicular to the plane containing the incoming and outgoing 
directions. Dipole radiation emits preferentially perpendicular to the direction of oscillation, so 
the light we observe will be polarized. To get the total effect we need to integrate over the light 
intensity in all the directions. If the radiation incident on the electron is isotropic there will 
be no net polarization generated after the scattering, simply by symmetry. Dipole distribution 
also does not generate polarization. For example, consider the case where the intensity of the 
incident radiation is higher from the top and lower from th e bottom, with the average intensity 
incident from the sides. Observing from above the plane scattered light from photons incident 



from either the bottom or top of the plane will be polarized in the horizontal direction, while 
that coming from the sides will be polarized in the vertical direction. Adding up all the compo- 
nents we find that horizontal and vertical components are equal, implying no polarization. We 
therefore need a quadrupolar pattern of intensity, in which case there is an excess of horizontal 
polarization from both top and bottom, which does not cancel the smaller vertical polarization 
from the sides. 

Since we need Thomson scattering it is clear that polarization cannot be generated after 
recombination of electrons and protons into hydrogen, when the universe becomes transparent 
for the photons. We still need to show however that there is sufficient quadrupole moment prior 
to last scattering, because only then will polarization be generated. In the electron rest frame 
one sees photons coming from distances of the order of a mean free path Xp. The photons that 
are scattered off a given electron come from places where the fluid has velocity v and because of 
the tight coupling between photons and electrons the photon distribution function has a dipole 
term Ti = n • v. To generate quadrupole moment one requires a gradient in the velocity field 
across the mean free path, T2 = Xpri^n^ diVj in the rest frame of the electron. But since the mean 
free path was so short prior to recombination the induced quadrupole was very small, until the 
last moment during recombination when mean free path started to grow. This implies that the 
mean free path is of the order of the duration of recombination and a detailed analysis gives 1151 

(Q + i[/)(n) « 0.17Anm*m^9it;j|,., (6) 

where Ar* is the width of the last scattering surface, happening at r*. The appearance of 
m = ei + ie2 in equation ^ assures that {Q + iU) transforms correctly under rotations of (ei, 62) 
(equation [1^. 

Figure y shows the temperature, E polarization and T — E cross-correlation power spec- 
tra. The oscillations in the spectra can be understood with the tight coupling approximation 
described above. At the time of recombination waves of different wavevectors k are at different 
phases in their oscillations. The phase only depends on the amplitude of the wavevector, so all 
the modes with the same amplitude of wavevector have the same phase. For particular values 
of this amplitude all the corresponding modes vanish. A given physical scale is related to the 
corresponding angular scale via the angular diameter distance, but the correspondence is not 
perfect because of 3-d to 2-d projection. If it were perfect the power spectrum of polarization 
would vanish at some values, but in reality it does not (figure ^). Nevertheless, the oscillations 
in polarization are much more pronounced than in temperature, because the latter receives con- 
tributions both from density and velocity and the two are out of phase with each other, so that 
even if one vanishes for some value of k the other does not. 

On large scales polarization is strongly suppressed. Correlations over large angles can only 
be created by the long wavelength perturbations, but these cannot produce a large polarization 
because of the tight coupling between photons and electrons prior to recombination: only wave- 
lengths that are small enough to produce anisotropics over the mean free path of the photons will 
give rise to a significant quadrupole in the temperature distribution, and thus to polarization. 
On small scales (/ > 1000) both polarization and temperature anisotropics are suppressed by 
photon diffusion. Here the wavelength of the perturbation becomes comparable to the mean free 
path of the photons prior to recombination. Photons can diffuse out of density peaks smaller 
than a diffusion length, thereby erasing the anisotropics. Note also that there is an extra gradient 
in polarization relative to temperature (equation which explains why polarization amplitude 
peaks at a smaller angular scale. The fact that the polarization field has relatively more small 
scale power is particularly evident when we compare the the T and E fields in plates 1 and 3. 
Finally, temperature polarization cross-correlation oscillates just like the other two spectra, but 
can be either positive or negative (figure |3|). 

Figure § shows the correlation functions in real space. The spectrum has been smoothed 
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Figure 3: Power spectra for COBE normalized SCDM. The E and the T — E spectra have been rescaled for 

convenience. 

with a Ofwhm = 0.2° gaussian, similar to the MAP beam. The correlation functions are defined 
in their natural coordinate system. In this frame there is no cross-correlation between Ur and T. 
An interesting point is that both polarization auto-correlation functions are negative for some 
angular range, which does not happen for the temperature. To interpret the cross correlation 
we can consider the polarization pattern around a hot spot (T > 0). The cross correlation 
starts positive, implying a radial pattern of polarization. Not all the polarization is correlated 
with the temperature so it is difficult to see this trend in plate 3. In plate 4 we only plot the 
correlated part of the polarization and here one can see more clearly that polarization "vectors" 
are preferentially radial around hot spots. 

Moving to larger angular scales the cross correlation changes sign. When it is negative the 
pattern becomes tangential. For large separations the polarization around a hot (cold) spot 
is tangential (radial). A point worth noting is that the cross correlation vanishes at = 0, in 
contrast to what happens for the other correlation functions. This follows simply from symmetry, 
since there is no preferred direction that polarization should take. Only when we consider 
two points separated by some distance is the symmetry broken. The vector joining the two 
points becomes the privileged direction and the polarization can be preferentially parallel or 
perpendicular to this direction. 

4 Conclusions 

Polarization has a rich geometrical structure, which can be simply understood using a real space 
construction of scalar and pseudoscalar fields. These can be constructed as integrals over Stokes 
Q and U parameters. Finite extent versions exist which allow one to search for B polarization 
without measuring the whole sky. 




Figure 4: Correlation functions in real space for COBE normalized standard CDM model. The spectra have been 
smoothed with a Bfmhm ~ 0.2° corresponding to the beam size of MAP. 



Generation of polarization requires both Thomson scattering and significant quadrupole 
moment of photon distribution in electron rest frame. We give simple arguments why this is 
so and discuss predictions for some realistic models. Realistic attempts to extract polarization 
information have to include complications such as instrument noise and foregrounds and are 
given elsewhere in these proceedingslHI. 
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